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■ Abstract 

^■i We report some improved wave function for mesons taking linear confinement term in 

1 standard QCD potential as parent and Coulombic term as perturbation while applying 

quantum mechanical perturbation technique in solving Schrodinger equation with such a 
potential. We find that Airy's infinite series appears in the wave-function of the mesons. 
We report our calculations on the Isgur-Wise function and its derivatives for heavy-light 
mesons, within this framework. 
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1 Introduction: 



psj . Understanding the physics of hadrons, specially the heavy-light meson, has been the focus 
of interest in recent past. For heavy-light meson, the heavy-quark mass is much greater than 
h-QCD, the QCD scale parameter. The four velocity of the heavy-quark is essentially the same 
as the four velocity of the meson and the heavy-quark, in the rest frame, appears as a static 
colour source. This gives rise to the consequence of spin-flavour symmetry in HQET and form 
factor suppression [1] . The overlap between the state of the light degrees of freedom in presence 
of a colour source moving with a velocity v , and that in the presence of a colour source moving 
with velocity v', is the Isgur-Wise function(IWF)[2]. This is the fundamental quantity in QCD, 
^ ■ which can be determined non-perturbatively. 

The motivation behind the calculation of IWF is to be able to extract CKM matrix elements 
V c b [3] from heavy-quark decays. Of particular interest are the semi-leptonic decays [4] like 
B — > DlV and B — > D*lT> , from which it is possible to extract Vcb- 

The main part of the IWF is the wave function of the meson [5] and some kinematic factor 
which depends upon four velocities of heavy-light mesons before and after recoil. Thus, the 
construction of wave function of meson is important for the detailed study of IWF and its 
derivatives like slope (charge radii of meson) and curvature (convexity parameter). 
For obtaining the meson wave function in potential model approach, the choice of potential 
is the first step. There are several potential models for quark-antiquark bound states like 
the Martin Potential [6], Cornell Potential [7], Richardson potential [8], Logarithmic potential [9]. 
The basic condition followed in enunciating these potentials are their flavour independence and 
existence of linear confinement. Of all the power law potentials, the Cornell potential is be- 
lieved to be the more realistic phenomenological potential model for mesons as it is based upon 
extrapolation of two kinds of asymptotic - ultraviolet at short distance (Coulombic term) and 
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infrared at large distance (linear term). It thus includes both the QCD concepts of 'confine- 
ment' and 'asymptotic freedom'. 

Starting from this Cornell potential of linear plus Coulombic form [10] and using two-body 
Schrodinger equation, several successful attempts have been reported in recent years in extract- 
ing the wave function of heavy-light mesons . This has been done both with linear confinement 
term in potential as perturbation [11] and Coulombic term as perturbation [14]. 
The case of Coulombic part in potential as parent in perturbation technique being simpler, 
several successful attempts have been made earlier [10-13] with such a choice, to extract meson 
wave function and make subsequent studies of meson properties. However, the origin of qq 
potential and also lattice QCD support the fact that the linear confinement term is the main 
contributing part to the potential. This inspires us to consider linear part in potential as parent 
and Coulombic part as perturbation. 

The linear part in potential as parent leads to Airy's infinite polynomial function [16] in unper- 
turbed wave function. This brings some constraint in the perturbation technique in deducing 
the perturbed wave function with Coulombic part as perturbation. Some approximated form 
of the perturbed wave function has been reported recently [14] considering Airy function up to 
0(r 3 ) only, which obviously has its limitations. 

In this paper, we report improvement of the previous wave function [14], and construct the 
meson wave function (with Coulombic part in potential as perturbation) containing complete 
Airy's infinite series, following Dalgarno's perturbation method [17]. 

Analysis of IWF with such a wave function containing Airy's infinite series, has been found to 
bring divergences. Very recently, we have overcome the problem of integrability of the otherwise 
divergent infinite Airy's series in IWF, by successfully introducing some reasonable cut-off to 
the infinite upper limit of integration [15]. 

In this paper also, we follow the same approach with our improved wave function and consider 
some cut-off to the infinite upper limit of integration while studying IWF and its derivatives. 
We calculate some reasonable range of this cut-off value which best suits our results, when 
compared with theoretical and experimental results for derivatives of IWF. We also study the 
sensitivity of the different orders of polynomial approximation of the infinite Airy's function in 
our calculation. 

In section 2 we discuss the necessary formalism while section 3 contains the calculations and 
results. Summary and concluding remarks are stated in Section 4. 



2 Formalism: 

2.1 Potential Model: 

The QCD inspired potential model of linear plus Coulombic type, has been reported in ref.[10- 
15]. For our studies of meson properties, we choose the widely accepted Cornell potential which 
has the form: 

V(r) = -C F — + br + c (1) 

r 

Cf is the colour factor, which is given by : 
N c is the colour quantum number; for N c = 3, we have Cp — I- 
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Here we take br as parent so that our unperturbed Hamiltonian [18] is: 

= + * (3) 

with 

H ' = ~t + C < 4 > 

as perturbation. Here \i is the reduced mass of the meson, which is given by: 

^=^^ (5) 
m q + mQ 

In the infinite quark mass limit (mg — > oo), /z ~ m 9 . We take b = 0.183 GeV 2 from charmo- 
nium spectroscopy [19] and c = 1 Gel / [20]. Under this consideration, the two body Schrodinger 
equation [21,22] for the Hamiltonian H = H + H' is: 

H\V >= (H + H')\V >= E\V > (6) 

2.2 Wave Function: 

To find the unperturbed wave function corresponding to H Q , we use the two-body radial 
Schrodinger equation (with h — 1): 

where we have put V(r) = br. 

For / = state, introducing u(r) = rR (r) and dimensionless variable g, we use the following 
formalism from ref. [23]: 

e =W*r-(fy/*E (8) 
The above equation (7) then reduces to: 

cPu , . 

^-m = o (9) 

The solution of this second order homogeneous differential equation [24] contains linear combi- 
nation of two types of Airy's functions Ai[r] and Bi[r]. From the nature of the Airy's function 
[25] we find that as r — > oo, Ai[r] — > and Bi[r] — > oo. So, it is reasonable to reject the Bi[r] 
part of the solution. The radial wave function thus has the form: 

u(r) = NAi[(2^b) 1 / 3 (r-j)] (10) 

where is our normalization constant which has the dimension of GeV 1 / 2 . The boundary 
condition u(0) = [26] gives us the unperturbed energy for ground state [23]: 

w° = E = -(^ 3 eo (ii) 

where £ is the zero of the Ai[r], such that Ai[^ ]=0 [24] and it has the explicit form: 

= _[5Mipi) r /3 (12) 
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In our case we consider ground state wave function only (n = 1). 
The unperturbed wave function for ground state is thus obtained as: 

tt°(r) = ^Ai^b) 1 ^ + g ] (13) 
*o (r) = ?L Ai[gir + ^ = -^—Ai\q\ (14) 
where we have taken g 1 = (2 fib) 1//3 and £>(r) = g±r + £> - 

The first order perturbed Eigen function \&'(r) and Eigen energy W can be calculated using 
the following relation: 



Hq$' + #'* = W°W + W'V° (15) 



We find, 



CO 

W = I r 2 H' | ^° | 2 dr (16) 



o 



Now taking H' = — f with B = we obtain from equation (15), 

+ H ] + br ~ EW{r) = [ 7 + w ' mr) (17) 

While applying Dalgarno's method of perturbation [27], we introduce the identity: 

Ai!{g) = = Z(g)Ai(g) (18) 

so that 

Ai"(g) = Z 2 (g)Ai(g) + Z'(g)Ai(g) (19) 

Here, Z(g) is a function of r, which can be calculated from the explicit polynomial expansion 
of Airy's function [Appendix-B]. 

We obtain, after employing Dalgarno's perturbation method [Appendix- A] , 
N' 

tf'(r) = — L4i(r)r + A 2 (r)r 2 + A 3 (r)r 3 + A 4 (r)r 4 + A 5 (r)r 5 + ]Ai[g 1 r + g ] (20) 



Here, Ai(r), A 2 (r), A 3 (r) etc are obtained as : 

1\iB 
2g 1 k l + gjk 2 
2fiW 
2 + 4giA;i + gjk 2 



^) = - o I 31 (22) 



12 + 8ft*! + effca ( j 

A5(r) " 20 + lOftfcx + rffe (25) 
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where we take k\{r) and k 2 (r) to be some function of r, as given below : 



Z{q) = ^ (26) 
k 2 (r 



Z\g) + Z\g) = (27) 



r 



Thus, we obtain the total wave function , as: 



N' 

¥ otal (r) = — [1 + A x (r)r + A 2 (r)r 2 + A 3 (r)r 3 + ]Ai[ Ql r + g ] (28) 

Here N' is the normalization constant of total wave function which also has the dimension of 
GeV 1 ! 2 . Considering relativistic effect on the wave function, the total relativistic wave function 
is given by [21]: 

N' r 

= —[1 + M{r)r + A 2 {r)r 2 + A 3 (r)r 3 + \Ai[ Ql r + g }(— )~ £ (29) 

r a 

Here, 



"^i^rw^^'-y 1 -^^ 1 -^^ (30) 

The equations (20), (28) when compared with equations (16), (21) of ref[15], reveal the improve- 
ment of the present formalism over the previous one. Earlier, the perturbed wave-function was 
limited up to Airy function of order r 3 , whereas in the present formalism, the perturbed wave 
function is the product of complete Airy's polynomial series and another infinite series. 



2.3 Isgur-Wise Function: 

Under Heavy Quark Symmetry ( HQS ), the strong interactions of the heavy quarks are in- 
dependent of its spin and mass [28] and all the form factors are completely determined, at all 
momentum transfers, in terms of only one elastic form factor function, the universal Isgur-Wise 
function £(v,v'). £,(v,v') depends only upon the four velocities v v and v v > of heavy particle 
before and after decay. This £(v,v') is normalized at zero recoil [29]. If y = v u .v u / , then, for 
zero recoil (y — 1), = 1. With increasing recoil y grows. In explicit form IWF [30] can be 
expressed as: 

^y) = l-p\y-l) + C(y-lf + --- (31) 

p 2 is the slope parameter at y — 1, given by : 

P 2 = (32) 



Sy 



p is known as the charge radius. 



C is the convexity parameter given by : 



C-—\ y=l (33) 



The calculation of this IWF is non-perturbative in principle and its depends upon the meson 
wave function and some kinematic factor, as given below : 

poo 

£(y) = / Anr 2 \^(r)\ 2 cos(pr)dr (34) 
Jo 
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where cos(pr) = 1 — + ^ — h with p 2 = 2fi 2 (y — 1). Taking cos(pr) up to 0(r 4 ) we 

get from equation (34), 



£(y)= A-nr 2 \m{r)\ 2 dr-[ATi^ r 4 \f(r)\ 2 dr](y - 1) + [-tt/x 4 / r 6 |*(r)| 2 dr](y - 1) 2 (35) 
Jo Jo 3 Jo 

Equations (31) and (35) give us : 

/•oo 

p 2 =4tt^ 2 / r 4 |^(r)| 2 dr (36) 



o 







2 Z" 00 

C=-nfj, A r 6 |^(r)| 2 dr and (37) 

3 Jo 

47rr 2 \y(r)\ 2 dr = 1 (38) 
Equation (38) gives the normalization constants N and N' for \l/ (r) and ^lf t {r) . 

3 Calculation and result: 

From equation (29), we find that the total wave function contains two infinite series - one 
a power series in r with coefficients Ai(r), ^(r) etc and the other an infinite Airy's series 
(equation(B.l) in appendix-B). As the infinite limit of integration in IWF brings divergence 
of both the series, we are compelled to consider some reasonable cut-off tq to this upper limit 
of integration. In principle, this cut-off ro should be greater than the size of hadrons. In our 
calculations, we consider this r to be greater than rrih being the mass of hadron. This 
always keeps r to be much greater than size of meson concerned. Also, consideration of such 
cut-off to upper limit of integrations will not sacrifice the nature and value of IWF and its 
derivatives, because, Airy's function falls very sharply (almost exponentially) and almost dies 
out with increasing r- value [31]. In fact, the Airy's function value becomes negligibly small for 
r > 4 ( AiryAi[A] = 0.000952). 

Also, the graph of normalization constant (N' )versus the cut-off to upper limit ( r ) shows that 
N' value decreases with increase in r . For D meson, as a representative case, such variation is 
shown in fig 1. 

ki and k 2 in equations (21-25) are calculated considering lowest Airy polynomial order, as 
[Appendix-B] : 

h(r) = l + ^ (39) 

k 2 (r) = ^ (40) 

with k = ^^o (41) 
biQi 

Here, ai and b\ are constants in Airy's infinite polynomial with a\ = 3 2/3^ 2 /3) = 0-3550281 and 
h = 3l/3l ! (1/3) = 0.2588194. 
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Table 1: p 2 and C for D meson ( p = 0.2761) 



r (in GeV r ) 


N' p 2 C 


5.0 


0.6616 0.6821 0.1138 


5.5 


0.5249 1.2249 0.3571 


6.0 


0.3655 1.9099 0.7342 


6.5 


0.2566 2.463 1.1247 



Table 2: p 2 and C for B meson ( p = 0.318) 



r (in GeV x ) 


N' p 2 C 


4.0 


0.7277 0.6699 0.0979 


4.5 


0.7184 0.7030 0.1118 


5.0 


0.6069 1.1758 0.3451 


5.5 


0.4172 2.0523 0.8722 


6.0 


0.2792 2.7854 1.4441 



Ai(r) and A 2 (r) etc as in equations(21-25) then have more explicit forms like: 



-2pBr 2 
2q x t 2 + 2g 1 kr + g\k 2 



Mr) = — T,L («) 



-2pW'r 2 

Mr) = (2 + 4 Ql )r 2 + 4 Ql kr + g 2 k 2 ^ (43) 



With these more explicit functions of r in the meson wave function as in equation (29), we have 
explored £(y) and its derivatives for Airy polynomial order up to r 10 , taking different cut-off 
values ranging from r = 3 GeV^ 1 to r = 7 GeV^ 1 , for D, D s , B, B s mesons taking a s = 0.22 
[32] . The results are shown in Tables 1 to 4. 

We find that, for the given Airy order, with increase in cut-off value, p 2 and C values 
increase steadily. However, the results show closer resemblance to recent results [33-41] ( as 
shown in Table 5) up to a specific cut-off value r for different mesons. For such specific range 
and order, our results show improvement over the result of ref [8]. At cut-off value higher than 
r ~ 6 GeV' 1 , the results jump to higher values than our expectations ( Figure 2(a) and 2(b)). 

Regarding sensitivity of the order of polynomial in infinite Airy's function, the results for 



Table 3: p 2 and C for D s meson ( p = 0.368) 



ro (in GeV r ) 


N' p 2 C 


4.4 


0.7648 0.7968 0.1388 


5.5 


0.7014 1.0707 0.2882 


5.0 


0.4947 2.1054 0.9665 


5.5 


0.3145 3.1059 1.8176 


6.0 


0.2121 3.8782 2.6875 
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Table 4: p 2 and C for B s meson ( p, = 0.4401) 



r (m GeV x ) 


TV' p 2 C 


3.0 


0.8511 0.8792 0.1651 


3.5 


0.8087 0.9809 0.2101 


4.0 


0.7855 1.0924 0.2805 


4.5 


0.5825 2.2312 1.1423 


5.0 


0.3512 3.6471 2.5282 



Table 5: Results of slope and curvature of £(y) in different models and collaborations. 



Model / collaboration 


Value of slope 


Value of curvature 


Ref [14] 


0.7936 


0.0008 


Le Youanc et al [33] 


> 0.75 


> 0.47 


Skryme Model [34] 


1.3 


0.85 


Neubert [35] 


0.82±0.09 




UK QCD Collab. [36] 


0.83 




CLEO [37,38] 


1.67 




BELLE [39] 


1.35 




HFAG [40] 


1.17 ±0.05 




Huang [41] 


1.35 ±0.12 





p 2 and C do not differ much upon variation of order of polynomial in Airy's function from r 4 
to r 10 . The results of such study in the case of D meson is shown in Table 6. 

In all calculations, we have taken care that the boundary condition of IWF ( £(1) = 1 ) is 
maintained throughout. The variation of £(?/) with y for different cut-off value ro ,with Airy 
order r 10 , for different mesons is shown in Figures 3(a) to 3(d). 

4 Conclusion and remarks: 

In this work, we have deduced the wave function for mesons involving the complete Airy's 
polynomial series, without introducing any approximation. The wave function constructed 

Table 6: Sensitivity of the polynomial order of Airy's function( for D meson with ro = 5.48) 



Airy Order 


P 2 


C 


r 4 


1.6338 


0.3345 


r 6 


1.1657 


0.3346 


r 7 


1.1666 


0.3347 




1.1665 


0.33467 


r lu 


1.16644 


0.334674 


With complete Airy series 
(After numerical integration) 


1.16644 


0.334674 



8 



Table 7: Values of p 2 and C with asymptotic form of Airy's function. 



r (in GeV- 1 ) 


p 2 (asymptotic ) 


C (asymptotic) 


3 


8.4 x 


1.2 x 10- 6 


4 


2.6 x 10~ 7 


6.0 x 10~ 8 


5 


4.6 x 10" 9 


1.6 x 10" 9 


6 


5.027 x 10" 11 


2.464 x 10" 11 


7 


3.56 x 10" 13 


2.345 x 10" 13 


8 


1.695 x 10" 15 


7.028 x 10" 15 



here is certainly an improvement over that in references [14,15]- in which the perturbed wave 
function W(r) has been calculated using Airy order r 3 . Our perturbed wave function(and hence 
total wave function) contains complete Airy's series. 

In our earlier work [15], we have shown that cutting off the upper limit of integrations in £(y) 
and its derivatives to some reasonable point does not upset the result, rather it almost conforms 
to the experimental expectations. In this work, we have followed similar approach in the study 
of IWF. Such study also, in turn, confirms the compatibility of our formalism. 
Also, for each value of cut-off r , we have considered the asymptotic form of the Airy's function 
taking limits of integration from r to oo . 

exD f_2 3/2\ 

- (44) 

With this asymptotic form we have also calculated the derivatives of £(y). Such analysis shows 
that, taking this asymptotic form of Airy's function very small values of p 2 and C are ob- 
tained(Table 7). This, thus, also confirms that the margin of error, by considering cut-off to 
infinite upper limit of integration, is negligible. 

The figures 3(a) to 3(d) confirm the fact that boundary condition for zero recoil (£(1) = 1) 
is maintained all through, with given polynomial orders of Airy's function and for different 
cut-off values. 

Regarding the results of ref [14] , which is for Airy order r 3 , we find that results match with 
our cut-off value r = 5.095 GeV~ x for D meson. 

Further, we study the compatibility of our potential model with the recent results of Heavy 
Flavour Averaging Group(HFAG) [40]. Taking the result of p 2 , we fix the value of cut-off for 
different orders of polynomial in Airy function (Table 8). This indicates that, the range of 
cut-off value r = 4 GeV' 1 to r = 6 GeV' 1 matches the expectations of ref [40]. 



W conclude by making the following comments on our present approach. 

• We opt for the quantum mechanical perturbation technique in deducing the meson wave 
function, due to the constraint in getting exact analytic solution of Schrodinger equation 
involving linear plus Coulombic type potential. 

• It is well known that method of numerical solution of Schrodinger equation [42] or Lat- 
tice QCD [43] gives more accurate results compared to relatively crude potential model 
approach following perturbation technique. However, the potential model approach is 
believed to give more physical insight into the problem under consideration. 
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Table 8: Value of cut-off r fixing p 2 = 1.17 from ref[40] 



meson 


r (in GeV x ) 


D 


5.48 


B 


4.998 


D s 


4.54 


B s 


4.05 



• Lastly, in our analysis of IWF, the parameters k\(r) and ^(r) and hence Ai(r),A 2 (r) 
etc are calculated considering lowest order in Airy polynomial. Considering higher Airy 
polynomial order in such calculations may bring further refinement of the result. 
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A Appendix 

From equation (17), in terms of radial wave function R{r) we get: 

[-^ + ^) + ^-^« = [7 + ^M (A.i) 

Let, 

R(r) = ^F(r)Ai[g\ = ^-F^Ai^r + g ] (A.2) 

So that, 

^ = -±-F(r)Ai[g] + -F'(r)Ai[g] + ^F(r)Ai'[g] (A.3) 
dr r l r r 

and 

^ = ^F(r)Ai[g]-^FXr)Ai[g]-^F(^ 

(A.4) 

Now we introduce the identity : 

Aif[Q\ = = Z(g)Ai(g) (A.5) 

so that 

Ai"(g) = Z 2 (g)Ai(g) + Z'(g)Ai(g) (A.6) 
Then the equation (A.I) becomes : 

^F'\r) + ^F'(r)Z(g)+^F(r)Z'(r) + ^F(r)Z 2 (r)-2fi(br-E)^F(r) = -2/x(^+^) (A.7) 
so that, 

F"(r) + 2eiF'(r)Z( e ) + ^[^ 2 (^) + Z'(g)]F(r) - 2^(br - E)F(r) = - — - 2/iW (A.8) 

r 



Assuming 



= ^ (A.9) 



r 



^) + Z'(,) = M^ (A.10) 

We get, 

F"(r) + 2 ft F'(r)^i + g\F{r)^- - 2fi{br - E)F(r) = - 2#' (A.ll) 

We take, 

F(r) = J2 A i rl ( A - 12 ) 

F'(r) = /^A^- 1 (A.13) 
i 

F"(r) = 1(1-1) Air 1 " 2 (A.14) 
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so that equation (A. 11) becomes, 



1(1-1) A l r l - 2 +2g 1 k 1 l ^ A^+gjh £ A,r'- 2 +2/iE ^ ^ A t r l+1 = ^-2^' 

z z / z z r 

(A.15) 

Equating powers of r 2 on both sides of equation (A.15) we get : 

g 2 1 k 2 A = (A. 16) 

This gives A = 0. 

Further, equating powers of r -1 , r° , r 1 , r 2 and r 3 of both sides of equation (A.15), we get: 

2g 1 k 1 A 1 + g\k 2 A x = -2fiB (A.17) 

2A 2 + A Ql hA 2 + g\k 2 A 2 = -2/jW (A.18) 

6A 3 + 6 gi hA 3 + £ 2 A; 2 A 3 + 2 f iEA 1 = (A.19) 

12A 4 + %Qik x A A + ^ 2 A; 2 ^4 + 2fiEA 2 - 2fiBA 1 = (A.20) 

20A 5 + IO^/ciA + £ 2 /M 5 + 2/x£A 3 - 2^M 2 = (A.21) 

From these equations (A.17-A.21), we extract the expressions of Ai(r), A 2 (r), A 3 (r), A 4 (r) and 
A 5 (r) as: 



2fiB 
2g 1 k 1 + g\k 2 
2fiW 

2 + 4eifci + ^fc 2 



^"■) -- ++pr 

*(0 = ~ n . SI (A.23) 



Mr) = - 'T? 2 , (A.24) 
6 + 6^i fci + £>ffc 2 

Mr) = -^ E t- 2 ""t (A-26) 



12 + 8glfcl + £ 2 /c 2 

A(r) = - 2 ^-2^ 2 

20 + 10^^ + ^^ V ; 

(A.27) 



F(r) = ^ Air 1 = Ai(r)r + A 2 (r)r 2 + A 3 (r)r 3 + A 4 (r)r 4 + A 5 (r)r 5 + (A.28) 

z 

Ultimately , we get the perturbed wave-function \&'(r) as: 
AT' 

*'(r) = — [Ai(r)r + A 2 (r)r 2 + A 3 (r)r 3 + A 4 (r)r 4 + A 5 (r)r 5 + ]Ai(^r + g ) (A.29) 
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B Appendix 

The Airy's infinite series as a function of g = Qir + g can be expressed as [24] : 

A[ ft r + go] = ai [l + + 18Q + 12960 + -] " 

uu . \ , {eir + QoY , (gir + gp) 7 (gir + g )10 

M(£ir + , ) + 12 + 5Q4 + 45 3 6 o + -1 (B-l) 

To find fci(r) and ^(r) , we take truncated Airy series up to lowest order, so that, we have 

fci(r) Ai'(^ir + ^ ) 



Therefore, 



and 



with 



Also, 



and 



so that, 

We thus obtain, 



Z(r) 



r Ai(gir + g ) 

oi - &i(£i r + go) 
_ Mi 
&i(gir + go) - ai 

= 1- _ gl ~ &lgp lj-i 

r bigx r 

-[l--]^ = -(l + -) 

/y* f* T* f* 



k 

fci(r) = 1 + - 
r 

*(0 = z + 4 



r r z 
bigi 



1 2A; A; 2 , 1 2A; A; 2 



B.2 
B.3 
B.4 
B.5 
B.6 



B.7 



B.8 



B.9 



^) = ^(l + - + ^) = ^ + -F + -I (B.10 



*M = -^-^ (b.ii 



^ = Z 2 (r) + ^(r) = ^ (B.12 



fca(r) = ^ (B.13 
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Figure 1: Variation of N' with ro 



pvsr 




7.5 
7.0 
6.5 
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4.5 
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2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 



(a) p 2 vs r 



(b) C vs r 



Figure 2: Variation of p 2 and C with r for different mesons 
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(a) £(y) vs y for D meson 



(b) £(y) vs y for B meson 



^(y) vs y for D meson 




£(y) vs y for B meson 




(c) vs y for D s meson (d) £(?/) vs y for i? s meson 

Figure 3: Variation of £(y) with y for different mesons 
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